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O . Abstract 

■ 

We estimate energy spectra and fluxes at the Earth's surface of the cosmic and 
Galactic neutrino backgrounds produced by thermonuclear reactions in stars. The 
extra-galactic component is obtained by combining the most recent estimates of the 
cosmic star formation history and the stellar initial mass function with accurate 
theoretical predictions of the neutrino yields all over the thermonuclear lifetime of 
. ^ ' stars of different masses. Models of the structure and evolution of the Milky Way 

^ ■ are used to derive maps of the expected flux generated by Galactic sources as a 

function of sky direction. The predicted neutrino backgrounds depend only slightly 
on model parameters. In the relevant 50 keV-10 MeV window, the total flux of 
cosmic neutrinos ranges between 20 and 65 cm -2 s _1 . Neutrinos reaching the Earth 
today have been typically emitted at redshift z ~ 2. Their energy spectrum peaks at 
E ~ 0.1 — 0.3 MeV. The energy and entropy densities of the cosmic background are 
negligible with respect to the thermal contribution of relic neutrinos originated in 
the early universe. In every sky direction, the cosmic background is outnumbered by 
the Galactic one, whose integrated flux amounts to 300-1000 cm~ 2 s . The emission 
from stars in the Galactic disk contributes more than 95 per cent of the signal. 
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1 Introduction 



Almost all our current knowledge of the universe derives from the detection of 
electro-magnetic quanta of different energies. However, the hot dense regions 
lying at the centre of most astrophysical sources are completely opaque to pho- 
tons, and their physical properties can only be inferred indirectly. Additional 
information on astrophysical objects is gathered by studying the properties 
of cosmic rays. Yet, these are charged particles and the presence of galactic 
magnetic fields makes nearly impossible to trace them back to their sources. 
Another limitation which applies to observations of both high-energy photons 
and charged particles is given by the so called Greisen-Zatsepin-Kuz'min effect, 
namely the creation of electron-positron pairs on the cosmic infrared radiation 
background (for photons) and of pions on the cosmic microwave background 
(for protons and heavier nuclei). In other words, the universe is optically thick 
to photons above ~ 10 TeV and to protons above ~ 10 EeV on a scale of few 
tens of Mpc. Neutrino astronomy, although still in its infancy, might provide 
a way to circumvent all these limitations. In fact, neutrinos are stable (so that 
they can cross long distances), weakly interacting (so that they are able to 
penetrate regions which are opaque to photons), and electrically neutral (so 
that their trajectories are not affected by magnetic fields). It is then reason- 
able to expect that neutrino telescopes would allow us to get spectacular new 
views of the universe. 

Regrettably, the same reasons which make neutrinos interesting astrophysical 
probes make their detection extremely difficult. At present, only the Sun and 
Supernova 1987A have been detected through their neutrino emission. This 
is not surprising, since, as like as in the electro-magnetic visible band, the 
Sun is expected to largely outshine any other astronomical source. Anyway, 
the detection of solar neutrinos has already led to important advances. By 
cross-checking the experimental results against the predictions of standard so- 
lar models it was realized that only a fraction (ranging between a third and 
two thirds depending on the particle energy) of the expected solar neutrino 
flux had been detected (see e.g. [1] and references therein). A possible inter- 
pretation is that part of the solar electron neutrinos transform into different 
particles. This indication, demonstrated by the Sudbury Neutrino Observatory 
(SNO) detection of non-electron neutrinos from the sun [2] and confirmed by 
the disappearance of reactor antineutrinos in the Kamioka Liquid Scintilla- 
tor Anti-Neutrino Detector (KamLAND) [3], requires modifying the minimal 
standard model of electro-weak interactions so as to allow for neutrino flavor 
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mixing and oscillations. 



An important problem for the development of neutrino astronomy is the de- 
termination of sky backgrounds as a function of energy. Many different sources 
might contribute, ranging from distant galaxies to Galactic supernovae. In this 
paper we estimate fluxes and energy spectra at the Earth's surface of the cos- 
mic and Galactic neutrino backgrounds generated by thermonuclear activity 
in stars. 

The nuclear fusion reactions which turn hydrogen into helium produce two 
electron neutrinos and release about 27 MeV of energy (including the neu- 
trino budget). The generation of a solar luminosity by hydrogen-core burning 
then implies the emission of about 2 x 10 38 neutrinos per second. In con- 
sequence, we expect the existence of a diffuse and nearly isotropic neutrino 
background originating from all stars in the universe that are and/or have 
been thermonuclearly active. At the same time, the fact that we live within 
the disk of a spiral galaxy should imprint characteristic features on the angular 
distribution of the neutrino flux. 

Despite the large effort aimed at estimating the amplitude and spectrum of the 
cosmic background generated by Type II supernovae [4,5,6] and Gamma-ray 
bursts [7], the current literature still lacks of accurate predictions concerning 
the cosmic and Galactic stellar neutrino backgrounds. A first attempt in this 
direction has been performed by Hartmann et al. [8] who, however, only con- 
sidered main-sequence stars and used semi-analytic stellar models to compute 
the neutrino emission rates. Accurate numerical determinations of the neu- 
trino yields over the whole H and He burning phases have been presented by 
Brocato et al. [9]. These results have been used by the same authors to esti- 
mate the present neutrino flux at the Earth's surface due to Galactic sources. 
Their Monte Carlo models give fluxes ranging from 3 to 24 cm -2 s -1 deg -1 
depending on Galactic longitude. A rough estimate of the cosmic background 
has been obtained by assuming that the present-day Milky Way is represen- 
tative of the whole galaxy population at all times. With this simple approach, 
Brocato et al. [9] found a background number density of n ~ 10~ 9 cm~ 3 and a 
flux c • n ~ 40 cm -2 s _1 (with c the speed of light in vacuum). As a matter of 
fact, the solar neutrino flux at Earth (~ 10 11 cm~ 2 s _1 ) is orders of magnitude 
larger than both the Galactic and the cosmic neutrino backgrounds. 

From the experimental point of view, it is presently impossible to discriminate 
the non-solar contribution from backgrounds in the detectors. For this reason, 
the cosmic and Galactic backgrounds will remain undetectable until techno- 
logical advancements will make possible the construction of highly-efficient 
directional neutrino detectors at energies ranging between 0.1 and 15 MeV. 

In this paper, the cosmic neutrino background is obtained by combining the 
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most recent estimates of the cosmic star formation history and the stellar ini- 
tial mass function with accurate theoretical predictions of the neutrino yields 
all over the thermonuclear lifetime of stars of different masses. We show that, 
even though this background is hardly of any cosmological relevance, its esti- 
mated amplitude, combined with observations of the 7-ray background, sets 
an upper limit to the radiative lifetime of neutrinos. By adopting state-of- 
the-art models of the structure and evolution of the Galaxy, we also derive 
maps of the expected flux generated by Galactic sources as a function of sky 
direction. Beyond contributing to the "noise" level for the detection of sin- 
gle neutrino sources, these backgrounds are themselves interesting since they 
carry information about the past history and structure of the Galaxy and our 
Hubble volume. Even though their detection is still beyond present-day tech- 
nical capabilities, detailed predictions of their properties may motivate and 
inspire future experimental efforts. 

The plan of the paper is as follows. In Section 2 we discuss the cosmic stellar 
neutrino background. After introducing the notation and summarizing the 
basic properties of thermonuclear neutrino sources, we present our results and 
discuss their cosmological relevance and the related uncertainties. In Section 
3 we estimate the amplitude of the Galactic neutrino background and derive 
maps of its variation on the sky. Our conclusions are listed in Section 4. 



2 Cosmic neutrino background 



2.1 Method 



Let us consider a homogeneous distribution of isotropic sources of neutrinos 
which fills the universe. Given their co moving neutrino emission rate e v (E,z) 
(particles per unit time, energy and comoving volume), the (mean) differential 
number flux of neutrinos (particles per unit time, surface, energy and solid 
angle) seen by an observer at redshift zq is 



I^\E) = 



dN v 



dtdSdEdQ 



c 
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+ z 
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dt 

dz 



dz , (1) 



(2) 



where 

— = H (1 + z) [fi (1 + zf + n k (1 + zf + n*] 1 ' 2 , 

with Qk = 1 — — ^a? fixes the relationship between cosmic time and redshift. 
In agreement with recent estimates (e.g. [10]), throughout this work we will 
assume that the present-day value of the matter density parameter fi = 0.3, 
the vacuum energy density parameter Q A = 0.7, and the Hubble parameter 
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H = 100/ikms Mpc with h = 0.65. It is convenient to separate the 
energy and redshift dependence in e v by writing, 

e v (E,z) = 52S®(E!)R$>(z) (3) 

i 

where the index % runs over the different thermonuclear channels for neutrino 
production. In this case, (E) is the spectral energy distribution of neutrinos 
at emission (normalized so that / S®(E) dE = l), while (z) gives the rate 
of neutrinos produced through the %— th channel. This can be estimated by 
combining a set of basic quantities. First, one needs to know when stars have 
been formed throughout the universe. This is parameterized by the cosmic star 
formation rate, ip(z), which gives the mass of stars formed per unit time and 
co moving volume as a function of redshift. Second, since stars with different 
masses produce different neutrino rates, it is important to know their relative 
abundances. This is described by the initial mass function (IMF), (f>(M), which 
gives the mass distribution of stars at birth. Finally, one needs to know the 
neutrino emission rate L$(M,t, Z) (particles per unit time) as a function of 
stellar mass, age, and chemical composition. In what follows, we will assume 
that the neutrino luminosity does not depend on the stellar metallicity, Z. This 
assumption is supported by numerical models showing that, independently of 
the stellar mass, varies by less than 10 per cent when the metallicity is 
changed from the solar value (Z Q ~ 0.02) to Z — 10~ 4 (see Ref. [9] and 
references therein). In summary, the intensity of the background observed at 
Earth due to the integrated effect of all the stellar sources in the universe is 



m E) = ^y r dz s<?>[(i + z)E\mz) (4) 

where the comoving emission rate of thermonuclear neutrinos is given by 

jm, n = jfcjjgjw A: ] dt> nz(f)} lw(m, t - f) 

iMZTdMM^M) ' { ' 

t(z) is the age of the universe at redshift z, and tu = max[0, t(z) — tuf e (M)] 
with tiif e (M) the lifetime of a star of mass M (intended as the time over which 
a star is producing energy through H-burning nuclear reactions). 



2.2 Stars as neutrino sources 



The net balance of hydrogen burning is the transformation of 4 protons and 
2 electrons into a 4 He nucleus. Conservation of lepton number requires the 
emission of 2 electron neutrinos. The energy yield of ~ 27 MeV is mostly 
released in the form of thermalized photons in the eV range (corresponding 
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to surface temperatures of 10 3 — 10 4 K). In other words, one neutrino leaves 
a star every 10 7 photons. The conversion of hydrogen into helium takes place 
through different channels which characterize the neutrino production rates 
and spectra for stars of different mass and age. Very-low mass stars (M < 
0.8M Q ) spend more than a Hubble time in their central hydrogen burning 
phase. They convert H into He mainly through the proton-proton I (ppl) 
chain. Neutrinos are emitted during the fusion of two protons into a nucleus of 
deuterium, p+p — * d + e + + u e , (pp neutrinos). Obviously, these stars had only 
the chance to produce pp neutrinos during their past lifetime. With increasing 
stellar mass, the ppll and ppIII chains become more efficient. In consequence, 
low-mass main-sequence (MS) stars emit a mixture of neutrinos generated by 
the pp reaction, the electron capture on 7 Be, and the 8 B decay. Note, however, 
that, in the subsequent H-shell burning phase, the energy production of these 
stars is dominated by the CNO cycle. In this case, also neutrinos coming from 
the f3 + decay of 13 N, 15 0, and 17 F contribute to the total flux. In more massive 
MS stars, the CNO cycle becomes progressively the dominant mechanism for 
energy generation, and for M > 1.5M Q CNO neutrinos dominate the MS 
stage. 

In this work, we consider all the nuclear reactions which produce electron 
neutrinos in the pp chains and the CNO cycle. This is done by using the 
stellar neutrino yields computed by Brocato et al. [9] as a function of stellar 
mass and age. These cover the whole H and He burning phases for stars with 
solar chemical composition and for nine selected values of the stellar mass 
ranging from 0.8 to 20 M Q . This set of evolutionary tracks covers the three 
stellar mass ranges characterized by different evolutionary behaviors related to 
the occurrence of electron degeneracy in He stellar cores (low-mass stars) or in 
C,0 stellar cores (intermediate-mass stars), or to the quiet ignition of nuclear 
reactions in non-degenerate cores (massive stars). Regrettably, no numerical 
model is available for M < 0.8 M . However, we found that for M < 1.7 M , 
the neutrino luminosity of pp neutrinos as a function of stellar mass is very 
well described by a power-law scaling relation oc M 4 4 . We use this function to 
extrapolate the emission rate of pp neutrinos down to 0.08 M (the minimum 
mass for hydrogen burning stars). 

While neutrino emission rates depend on the star structure, neutrino spec- 
tra are derived from the standard theory of electroweak interactions and are 
largely independent of stellar parameters [11] 1 . For pp and CNO neutrinos we 
use the energy distributions obtained in Refs. [13] and [14], respectively. The 
line shapes for the neutrinos emitted by the electron-capture process which 
causes the transition between the ground-state of 7 Be to either the ground- 
state or the excited state of 7 Li are taken from Ref . [12] . For the high-energy 



In principle, neutrino energy spectra are only slightly dependent on the temper- 
ature at which the nuclear reactions take place [12]. 
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tail of the neutrino spectrum due to 8 B decay we use the results presentend 
in Ref. [15]. 

In what follows, we will only consider electron neutrinos produced by ther- 
monuclear reactions burning H and He. In the advanced phases of stellar 
evolution, additional neutrinos and anti-neutrinos, of all three flavors, can be 
emitted directly at expenses of the thermal energy of the star (the so-called 
"cooling" neutrinos). They have energies between 10 and 60 KeV [16], and 
their expected emission rates are comparable to those of thermonuclear neu- 
trinos [9]. For this reason we will limit our discussion to energies larger than 
0.06 MeV. Oxygen and silicon burnings in massive stars (which respectively 
proceed at a temperature of about 160 and 270 keV) can easily produce ther- 
mal neutrinos up to the MeV scale. These processes are active for a short time 
(their duration is of the order of 1 yr for stars of 25 M ) and liberate a lot less 
binding energy per nucleon than hydrogen burning, so that the total energy 
release from such phases is significantly smaller. However, at variance with 
hydrogen burning, almost all of the energy released from these phases goes 
into neutrinos. Also, Type la supernovae are powerful neutrino sources in the 
sub-MeV range while, above a few MeV, core-collapse neutrinos are expected 
to dominate the counts. We defer the analysis of the background contribution 
of these sources to future work. 
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Fig. 1. Cosmic star formation history as a function of redshift. Different linestyles 
correspond to the functions listed in equation (6). 
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2.3 Star formation history 



In the past few years a number of authors have attempted to determine the 
star formation history of the universe from observational data. Although some 
details remain controversial, the evolution of the star formation rate (hereafter 
SFR) with redshift can be traced till to z ~ 6 using various diagnostics. It 
is widely accepted that for z < 1 the SFR is a rapidly increasing function of 
redshift, reflecting the fact that present-day galaxies have been assembled long 
ago. Corrections for dust reddening and extinction complicate the picture at 
higher redshifts. To summarize the present knowledge and the corresponding 
uncertainty, we use here four different parameterizations of the global star- 
formation rate per unit comoving volume: 



0.65 exp(2.93z) + 15 

, . . n ir 14exp[0.6(2-5.4)] ir llk „ 3 

M$ = 0-15 r I q n oU T7\T M © ^ M P C • 

5 + 9 exp|0. 933(2 — 5.4)J 

with 

(1 + ^)3/2 • 

The first three, obtained by Porciani & Madau [17], match present data from 
UV-continuum and H-a cosmic luminosity densities, and allow for different 
amounts of dust reddening, especially at high-z (see Ref. [17] for further de- 
tails). The fourth SFR is taken from Ref. [18] (see also [19]) and is based 
on theoretical speculation only. This is obtained from a set of cosmological 
simulations including gas hydrodynamics where an a priori recipe for star 
formation has been assumed. The whole analysis has been developed in the 
framework of a vacuum energy dominated cold dark matter model in which 
galaxies form hierarchically. At variance with the determinations based on ob- 
servational data, this SFR extends to very high redshifts but suffers from a 
number of theoretical assumptions. The resulting function has a maximum at 
z ~ 5.5 and declines roughly exponentially towards both high and low red- 
shifts. A semianalytic model to explain the origin of this functional form has 
been presented by Hernquist & Springel [20]. The SFRs given in equation (6) 
are plotted in Figure 1. The intensity of the cosmic neutrino background have 
been computed by integrating equation (4) out to ^ max = 10. 
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2.4 Initial Mass Function 



Basic reasoning suggests that the IMF should vary with the properties of the 
star-forming clouds. Nevertheless, although there is no fundamental reason for 
a universal mass distribution of stars, so far no firm evidence for a variable IMF 
exists (see e.g. Refs. [21,22,23] for different points of view regarding this issue). 
In what follows we will assume a time independent IMF, adopting different 
functional forms to take into account observational uncertainties. As a stan- 
dard model we use the IMF originally proposed by Salpeter: <p{M) oc M~ 2 - 35 . 
Even though this is supported by many studies, especially for M > 0.5 M & 
(see e.g. [24,25]), the behavior of the IMF at lower masses is still very uncer- 
tain (see e.g. [25,26,27]). This reflects the observational uncertainties in the 
determination of the star luminosity function at low luminosities. To represent 
a number of different observational results, we focus here on two additional ex- 
amples. From M-dwarf observations with the Hubble Space Telescope, Gould, 
Bahcall & Flynn [28] determined the following form for the low-mass end of 
the IMF (M < 1.6M ): 

log(0) = const. - 2.33 log(M/M ) - 1.82 [log(M/M )] 2 . (8) 

This function peaks at M ~ 0.23 M , and the slope at M — 1 M Q is —2.33, 
almost coincident with the Salpeter value of —2.35. We will call GBF IMF the 
distribution obtained by matching this equation for M < 1 M with a Salpeter 
IMF for higher masses. Finally, as a third option, we use a parameterization 
of the IMF obtained from star counts in different Galactic fields by Kroupa 
[25]. This is well fit by a broken power-law with a change in the exponent 
near 0.5 M/M. . namely, oc M~ ai , with a x = 1.3 for M/M Q < 0.5 and 
a 2 = 2.3 for M/M Q > 0.5. In all cases, we assume M min = 0.08 M and 
M mfac = 125 



2. 5 Results 

In this section, we present our main results, obtained by combining in equation 
(4) the different functions discussed above. In Figure 2, the expected cosmic 
background intensity, is plotted as a function of energy for different star 
formation histories and assuming a Salpeter IMF. The corresponding fluxes, 
integrated over all energies, are listed in Table l. 2 The amplitude of the 

2 In the literature, it has become customary to call flux the quantity cn = Airl v even 
though this does not correspond to any physical flux. In fact, for isotropically moving 
particles the net flux vector (n v) vanishes, while the (one-sided) flux through a fixed 
surface is given by F+ = J 2 ^I u (9,<f)) cosOdQ, which reduces to = ir I u in the 
isotropic case. 
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10- 2 0.1 1 10 

E (MeV) 

Fig. 2. Mean intensity of the cosmic neutrino background as a function of energy. 
Different linestyles correspond to different assumed SFRs. We only consider neu- 
trinos produced by thermonuclear reactions in stars. The shaded region highlights 
the energy range in which additional "cooling neutrinos" are expected to contribute 
significantly. 

background ranges within a factor of 2, with ip^ and ip\ giving the highest 
and the lowest fluxes, respectively. This is not surprising since ipi has a much 
lower amplitude than the other SFRs at z > 2, while ^3 keeps increasing with 
redshift. It is interesting to compare the contributions to the background due 
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Fig. 3. Contributions to the neutrino differential flux due to stars of different masses. 
A Salpeter IMF and 1^2 (2) are assumed. 
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Table 1. Cosmic neutrino background. 



Sr K 


Tl\ /TT? 

IMr 


cn 


n • 10 9 


u ■ 10 


\l lJ h ■ 10 


IP 

hi 




1 n7 /;,, 
S ■ 1U /Kb 






(cm^s- 1 ) 


(cm- 3 ) 


(MeVcm- 3 ) 




(MeV) 




(cm- 3 ) 


1 


Salp. 


20.2 


0.68 


2.73 


2.59 


0.40 


1.41 


0.61 


1 


Kroupa 


29.8 


0.99 


4.03 


3.82 


0.40 


1.45 


0.89 


1 


GBF 


35.9 


1.20 


4.86 


4.61 


0.41 


1.43 


1.07 


2 


Salp. 


26.0 


0.87 


3.29 


3.12 


0.38 


1.83 


0.77 


2 


Kroupa 


38.2 


1.27 


4.81 


4.57 


0.38 


1.88 


1.13 


2 


GBF 


46.1 


1.54 


5.85 


5.55 


0.38 


1.86 


1.36 


3 


Salp. 


36.7 


1.22 


4.67 


4.43 


0.38 


2.01 


1.09 


3 


Kroupa 


53.5 


1.79 


6.78 


6.43 


0.38 


2.06 


1.58 


3 


GBF 


65.0 


2.17 


8.32 


7.90 


0.38 


2.04 


1.92 


4 


Salp. 


21.2 


0.71 


2.56 


2.43 


0.36 


2.30 


0.63 


4 


Kroupa 


31.3 


1.04 


3.78 


3.59 


0.36 


2.36 


0.93 


4 


GBF 


37.6 


1.26 


4.57 


4.34 


0.36 


2.33 


1.11 



a This denotes the mean emission redshift for background neutrinos with E = 0.2 MeV. 



to low (M/Mq < 1.2), intermediate (1.2 < M/M Q < 8M ) and high mass 
stars (M/Mq > 8). This is done in Figure 3, where a Salpeter IMF and ^2(2) 
are assumed. The shape of the different curves can be explained as follows. 
Very-low mass stars, M/M Q < 0.9, produce only pp neutrinos with energies 
below 0.4 MeV. Masses in the range 0.9 < M/M Q < 1.2 pass through 
the ppll and ppIII chains, thus emitting 7 Be and 8 B neutrinos. The small 
contribution of 8 B neutrinos to the total flux is evident at higher energies. 
The peculiar feature at energies in the range 0.4 < (E/l MeV) < 2 is due 
to the reduced quantity of CNO neutrinos emitted by these stars during their 
H-burning shell phase. As discussed above, the stellar neutrino luminosity 
rapidly increases with mass and, in particular, the contribution of CNO, 8 B, 
and 7 Be neutrinos increases. Nevertheless, the steep slope of the IMF and the 
shorter main sequence lifetime diminish the contribution of massive stars to 
the stellar neutrino background. In summary, intermediate mass stars clearly 
dominate the background at all energies but for a narrow window around 1.5 
MeV where small mass stars give a slightly larger contribution. In Figure 4 we 
show how the spectrum of the background changes by varying the IMF while 
using the same SFR (in this case, ^ 2 ). The reader is referred to Table 1 for 
the corresponding fluxes. The net effect is a slight spectral distortion and a 
difference in the background intensity as large as a factor of 2. As expected, the 
amplitude of the background increases when the relative weight of low-mass 
stars is reduced. 
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Fig. 4. As in Fig. 1 but for different IMFs, assuming -02 for the SFR. 

It is interesting to study how the the probability density distribution for the 
emission redshift of neutrinos varies with their energy at detection. This is 
simply given by 



P(z,E) 



e v [E(l+z),z] 



1 



47T H [fi„(l + + fi fc (l + Zf + fi A ]V2 7 (o) (£7) 



(9) 
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Fig. 5. Probability distribution of the emission redshift for neutrinos with energy E 
in the cosmic background (solid line). Different linestyles denote the contribution 
of different generation processes, namely: pp fusion (dashed), 13 N decay (dotted), 
and 15 O decay (dot-dashed). Left and right panels refer to E = 0.2 MeV and E = 2 
MeV, respectively. 
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and it is shown in Figure 5 for two different neutrino energies. Note that 
neutrinos at the peak of the energy spectrum (E ~ 0.2 MeV) have been 
emitted within a wide redshift range. Three different generation mechanisms 
contribute: the nuclear fusion of 2 protons (which accounts for the 27 per cent 
of the neutrinos observed at E = 0.2 MeV) and the (3 + decays of 13 N (44 per 
cent) and 15 (29 per cent). Each of these processes corresponds to a different 
distribution of emission redshifts, which is determined by the joint effect of 
the SFR and the shape of the rest-frame emission spectrum. For instance, the 
narrow peak near z = is generated by stars with M ~ 1M© which started to 
produce energy through the CNO cycle only recently. We find that the mean 
emission redshift for photons at E — 0.2 MeV is z = 1.83. On the other hand, 
at higher energies, only the 8 B f3 + decay contributes. The emission-redshift 
distribution for neutrinos at E — 2 MeV is shown in the right panel of Figure 
5. In this case, we find z = 2.03. 

For E < 0.06 MeV, the contribution to the total neutrino fluxes of cooling 
neutrinos should eventually be taken into account. In fact, the estimated neu- 
trino fluxes at the Earth's surface due to plasma-neutrinos emitted by stars 
till the end of the Red Giant phase as well as during the last phase of cooling 
as white dwarfs is comparable to the pp-neutrinos flux [9]. 



2. 6 Cosmological relevance 



It is interesting to compare the energy and entropy densities of neutrinos 
generated by thermonuclear reactions in the universe with the corresponding 
quantities for the primordial neutrino background. The present-day number 
and energy densities of the cosmic neutrino background generated by stars are 



47T 



n 



An 



f ll 0) {E)dE and u = ^ J E (E) dE , (10) 

respectively. We define the entropy density for our non-thermal cosmic neu- 
trino background as 



s = 



4?r k h g s 
; /i 3 



jH In 



1 



f(E) 



1 



In [l-f(E)]\dE, 



where 



f(E) = h 3 c 2 



E 2 



6.36 x 10~ 41 (MeV 3 cm 2 s) 



E 2 



(12) 



fcb denotes the Boltzmann constant, g s is the number of helicity degrees of 
freedom, and h is the Planck's constant. This is derived from first principles 
in the Appendix and gives the standard result s = 4.202 k h n for a ther- 
mal background. Note that, since neutrinos in the cosmic background are far 
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from being degenerate (i.e. f(E) <C 1), equation (11) reduces, in practice, 
to s oc — / E 2 f(E) \n[f(E)] dE. Values for the number, energy and entropy 
densities of the background are listed in Table 1 for the different models we 
considered. Typically, n ~ lCT 9 cm~ 3 , it~5x lCT 10 MeVcm" 3 (correspond- 
ing to a contribution to the density parameter of VL V h 2 ~ 5 x 10~ 8 ) and 
s ~ 10 _7 /cbcm~ 3 . On the other hand, from the standard theory of neutrino 
decoupling in the early universe one expects a thermal distribution of relic par- 
ticles with T rel = 1.95 K. This contributes n re i = 56.5 cm -3 , fi rcl h 2 ~ 4.8 x 10~ 6 
(for massless neutrinos), 3 and s rc i = 203.5 cm -3 for each neutrino family 
(an equal contribution is given by the corresponding antiparticles) . We can 
thus draw the conclusion that thermonuclearly produced neutrinos hardly are 
of any cosmological relevance. In particular, their abundance and entropy are 
overwhelmed by the corresponding properties of relic neutrinos from the hot 
initial phases of the universe. However, since thermonuclear neutrinos from 
stars are characterized by much higher energies than the thermal relics from 
the big bang, the ratio of the energy densities of this two backgrounds is much 
smaller than the number density ratio. In both cases, the contribution to the 
curvature of space-time is anyway negligible. 



2. 7 Constraints on the neutrino radiative lifetime and electromagnetic form 
factors 

In the hypothesis of radiative neutrino decay, v — > V(_ + 7, cosmic neutrinos [y 
with mass m) decaying into some neutral fermion (possibly a lighter neutrino 
i/£ with mass me) over a Hubble time would be the source of a diffuse photon 
background. For a fixed energy E v of the heavy neutrino, the photon spectrum 
is continuum, with an end point at E 1 = E v {m 2 — m 2 t )/m 2 for ultrarelativistic 
neutrinos. Depending on the mass pattern, the photon energy can extend up 
to the MeV region. As derived in Appendix B, the background intensity of 
the resulting photons at z = is 

m = ^t^lim { 2 E) , (13) 

where r~ l is the radiative decay rate and 

f°° dz f°° 

° ( M (l + ,)Wl + ^(l + ^ + »A]^ / d/P ^ • 

a numerical coefficient of order unity. 



3 Assuming that the sum of neutrino masses amounts to, at least, 1/30 eV (as 
suggested by recent data) corresponds to Q re i h 2 > 3.5 x 10~ 4 . 
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Observations have shown the presence of a diffuse (high- latitude) photon back- 
ground in the MeV region that appears to be of extragalactic origin. Both a 
large number of unresolved point sources and truly diffuse processes might 
potentially contribute. However, the X-ray background (below 100 keV) is 
generally understood to arise primarily from the integrated emission of active 
galactic nuclei (such as Seyfert galaxies) with redshifts ranging up to 6 [29]. 
At slightly higher energies (E ~ 1 MeV), radioactivity in Type la supernovae 
is expected to be the most important contributor [30]. 

The amplitude of the diffuse 7-ray background in the MeV range is extremely 
hard to measure (both from the Earth and from space) due to the extraor- 
dinarily high atmospheric and instrumental background. In particular, the 
interaction of cosmic rays with the detector and nuclear line emission domi- 
nate the noise. Recent results from a number of satellite missions indicate the 
presence of an extragalactic component significantly lower than previous mea- 
surements. The Medium Energy Detectors of the A4 experiment onboard the 
High Energy Astrophysics Observatory 1 (HEAO-1) measured a background 
intensity j(E) = (2.62 ±0.05) x (E/0.1 MeV)~ 2 - 75±om cm" 2 s" 1 MeV" 1 si' 1 for 
0.1 < (E/l MeV) < 0.4 [32]. The Imaging Compton Telescope (COMPTEL) 
on the Compton Gamma Ray Observatory (CGRO) showed that the diffuse 
background at 0.8 < {E/l MeV) < 30 is consistent with power-law extrapola- 
tion from the HEAO result at lower energies [32], and provided a new fitting 
formula j(E) = (1.05±0.2) x 10~ 4 (E/5 MeV) _2 - 4±0 - 2 cm~ 2 s _1 MeV -1 sr _1 [33]. 
Finally, the Solar Maximum Mission 7-ray spectrometer (SMM/GRS) mea- 
sured a signal for 0.3 < (E/l MeV) < 8.5 which is ~ 20% lower than the 
COMPTEL determination [34]. 

If one attributes a fraction / of the diffuse background to cosmic neutrino 
radiative decay, one can derive a lower bound for the radiative decay lifetime, 
r 7 (the total lifetime is obtained by multiplication with the branching ratio 
for the radiative channel), 4 

> (2-6) x 1Ql2 g ^ (15) 



mc 2 f 

The radiative decay rate can be uniquely characterized by the magnetic and 
electric (dipole) transition moments, fj, vut and e uue , through (e.g. [35]) 



r 7 




(16) 

where fi 2 s = |/wj 2 + \^w e \ 2 and natural units (h — c—1) have been adopted 



4 This limit is obtained for photons of 0.1 MeV, corresponding to neutrinos at the 
peak of the cosmic spectrum. On the grounds of the recent experimental results on 
flavor mixing and oscillation, this lower bound is valid for any neutrino type. 
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for convenience. Assuming a mass hierarchy m ^> me and measuring the 
effective magnetic moment in units of the Bohr magneton, /i B , this corresponds 
to 



Note that, for mc 2 < leV (the favored mass range from current experiments), 
the radiative decay channel is vastly phase-space suppressed, and the phase- 
space factor m 3 in equation (16) strongly reduces the strength of the limit on 
/x e ff- Anyway, even by assuming / = 1, our upper bound is much more stringent 
than the limit [i e &/ fJ>B < 0.09 (mc 2 /l eV)~ 2 derived from the absence of decay 
photons of reactor V e fluxes [36], and from the solar neutrino flux, /x e ff//-*B < 
5 x 10~ 6 (mc 2 /l eV)~ 2 [37]. However, most of the hard X-ray background 
(2-10 keV) has been recently resolved into discrete sources [38]. The data 
compilation by Moretti et al. [39] shows that 90-95% of the background can 
be ascribed to discrete source emission. Accounting for scattered light, the 
resolved fraction might even equal 99% of the total [40]. Active galactic nuclei 
may similarly contribute even at 100 keV [41]. Assuming / ~ 0.01, makes our 
result comparable to the limit derived from the Supernova 1987A neutrino 
burst, /Xefr/z^B < 1.5xl0~ 8 (mc 2 /leV)~ 2 , (see Ref. [42] and references therein) . 



Our bound on /j, e g, however, is not competitive with other astrophysical de- 
terminations (see [44] for a comprehensive review) based, for instance, on 
plasmon decay in globular cluster stars (/les/l^B < 3 x 10~ 12 ) [45], the ob- 
servation of TeV 7-rays from Markarian 421 and 501 [46], helioseismology 
(pLes/fJ-B < 4 x 10~ 10 ) [44], and the primordial neutrino background (fu, e s/^B < 
10 -11 (mc 2 /l eV)~ 9 / 4 ) [47]. Laboratory limits on the diagonal and transition 
magnetic moments are typically obtained from measurements of the v e scat- 
tering cross section at low energy near a nuclear reactor. Current upper limits 
to the transition moments involving a particular flavor are: fi/ fiB < 1.0 x 10~ 10 
for v e [48], fi/fi B < 7.4 x 10~ 10 for [49], and /i//i B < 5.4 x 10~ 7 for v T [50]. 



3 Galactic neutrino background 

Will we ever be able to detect the cosmic neutrino background generated 
by stars directly? Beyond technological issues related to the availability of 
detectors with the required sensitivity, one needs to worry about the presence 
of other backgrounds. It is straightforward to think of neutrinos generated by 
thermonuclear reactions taking place within the Galaxy, since it is expected 
to have a similar spectrum to the cosmic one. In this section we compute the 

5 For Dirac neutrinos, much more restrictive limits on the magnetic moments can 
be derived from Supernova 1987A data [43]. 
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basic properties of this new background and its directional dependence on the 
sky. For analytical convenience, we will use here a different parameterization 
of the star formation rate, ip^(t'), denoting the mass of gas converted into 
stars per unit time in a given galatic component (thin disk, thick disk, halo 
and bulge) labelled by the index (j). The mean density of stars that are still 
thermonuclearly active today can then be expressed as 

!^dMM<t>{M) ' { ] 

where to — 14.5 Gyr is the present cosmic age in the adopted cosmological 
model. Therefore, the mean neutrino rate presently emitted by a star in the 
component (j) through the production channel (i) is 

Q(i,j) = iM-dM 0(M) J& dt> (f) Lg (M, t - t>) 
fM™dM<f>(M)ft t0 M dt>i;U)(t>) 

where, once again, tu = max[0, t(z) — tuf e (M)]. Note that, being a mean quan- 
tity, Gfrfi does not depend on the overall normalization of the star-formation 
rate but only on its temporal evolution. We account for the geometric struc- 
ture of the galaxy by introducing a set of functions, n^(v g ), which describe 
the present-day number density distributions of stars in the different galac- 
tic components as a function of the distance from the galactic centre. The 
differential neutrino flux detected on Earth is thus 

J<°> (r, E) = -L £ G™ S® (E) J dr n® (r + r ) , (20) 

where r denotes separations from Earth and r is our distance from the galactic 
centre. 

In the last decades, accurate star count surveys, gas-dynamical and stellar- 
kinematics studies, and the analysis of near-infrared brightness maps allowed 
us to identify the main structural features of the Milky Way and to formulate a 
number of standard density laws to describe its populations. Recently, a series 
of technological improvements made possible a more accurate determination 
of the corresponding free parameters. In this section, we will use these results 
to estimate the intensity and directional dependence of the Galactic neutrino 
background. In particular, for the outer galaxy, we will use the results from 
the high-latitude star counts obtained by the Sloan Digital Sky Survey [51] 
and by Siegel et al. [52] which homogeneously cover nearly 279 and 15 square 
degrees, respectively. Given the level of uncertainty, the design of the models 
is deliberately simple. Typically they account for two double-exponential disks 
(a thin and a thick one), a spheroidal halo and a central bulge. 6 Given the 



Note that the integration in equation (20) can be performed using standard nu- 
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degree of approximation of our calculations, we neglect the spiral structure 
in the disk. The corresponding best-fitting parameters are listed in Table 2. 
The results are strikingly similar but for the scale height of the thick disk, 
which, however, is poorly determined in both cases. Note that these smooth 
models are far from being perfect, and neglect many detailed features of the 
stellar distribution in the Galaxy. For instance, the outer halo might not be 
smooth but be composed by overlapping streams of stars (e.g. [53]). In general, 
residuals show that these models systematically overpredict stellar densities 
towards the galactic centre and underpredict them in the outer galaxy [52]. 

Because of the lack of low-latitude data, these models constrain the inner 
structure of the Galaxy only weakly. For instance, Siegel et al. [52] assume 
an r~ 3 density distribution for the bulge, while it has long been known that 
the central region of the Galaxy has a distinct stellar density profile that 
approximately follows an r~ 2 law (e.g. [54]). Recently, it has become widely 
accepted that our Galaxy is barred, as evidence accumulated over the last few 
years from a number of different observations. Along this line, a more accu- 
rate model for the light distribution in the inner Galaxy has been extracted 
by Binney, Gerhardt & Spergel [55] from the analysis of dust-corrected maps 
of the near-infrared brightness emission of the Galaxy (see e.g. [56]). The 
corresponding parameters are listed in Table 2. However, converting near-IR 
brightness maps into stellar densities its not an easy task since one has first 
to assume a mass-to-light ratio for the different components and then a stellar 
mass function. All these steps have a rather large error, and this leads to a 
large uncertainty in the final result. On the other hand, one could think about 
relating IR emission directly to either stellar density or, even, to neutrino lu- 
minosity. However, one needs to remember that a significant fraction of the 
IR brightness comes from dust and gas. Moreover, during the evolution out 
of the MS, the neutrino emission rate does not correlate anymore with the 
photon luminosity of a star, since the only source of thermo-nuclear neutrinos 
is a H burning shell, while the structure is partially supported by the triple 
a reactions. For simplicity, we assumed a constant conversion factor between 
infrared brightness and stellar number density in the different components. In 
the end, in order to bracket our predictions for the Galactic neutrino back- 
ground, we combine the data listed in Table 2 into two distinct models. The 
first one (Gl) is completely based on the results by Siegel et al. [52] and it 
is expected to give the highest neutrino flux. The second model (G2) is in- 
stead obtained by combining the disk+bulge fit by Binney et al. [55] with the 
halo+thick disk results by Chen et al. [51]. In both cases, whenever a given 
parameter is associated with a finite range of values in Table 2, we assumed the 
mean value of the quoted extremals. Both models are normalized by imposing 
that the number density of stars in the solar neighborhood belonging to the 



merical techniques with no need to use slowly converging Monte Carlo methods as 
in Ref. [9]. 
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Table 2. Structure parameters in the adopted Galaxy models 



Parameter 



Binney et al. (1997) 



Chen et al. (2001) Siegel ct al. (2002) 



Disk density law 

Disk scale length (pc) 

Disk scale height (pc) 

Disk local normalization (pc -3 ) 



Double exponential 
2500 
210 + 42 a 
0.1 



Double exponential Double exponential 
2250 2000-2500 
330 350 
0.1 0.1 



Thick-disk density law 
Thick-disk scale length (pc) 
Thick-disk scale height (pc) 
Thick-disk/disk local 
normalization (%) 



Double exponential Double exponential 



3500 
580-750 

6.5-13 



3000-4000 
900-1200 

6-10 



Halo density law 
Halo power-law index 
Halo axial ratios (b/a, c/a) 
Halo/disk local normalization (%) 



Power law 
2.5 
(1,0.55) 
0.125 



Power law 

2.75 
(1,0.5-0.7) 

0.15 



Bar/Bulge density law 
Bar/Bulge axial ratios (b/a, c/a) 
Bar/Bulge power-law index 
Bar/Bulge scale length (pc) 
Bar/Bulge core size (pc) 
Angle between the Bar/Bulge 

major axis and the Sun-centre line 
Bulge/disk local normalization (%) 



Exp-truncatcd power-law 
(0.5,0.6) 
1.8 
1900 
100 
20° 

6.1 x 10" 13 



Power- law b 
(1,0.5) 
3 



0.02 



Solar distance 

from the Galactic centre (pc) 
Solar distance 

from the Galactic mid-plane (pc) 



8000 



14 



8000 



27 



8000 
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a In this model the vertical structure of the disk is given by the linear combination of two exponential laws with 
different scale heights. See Binney ct al. (1997) for details. 

b The bulge structure has been fixed a priori in this model and it is not determined by the data 

thin disk is 0.1 pc~ 3 [57]. Moreover, in order to avoid unphysical divergences, 
we adopt a constant density distribution within the inner 20 pc of the bulge 
(in Gl) and within the inner 240 pc of the halo (in both Gl and G2). 

No compelling constraints currently exist on the star formation histories of the 
various galactic components. We adopt an overall picture based on the studies 
by Binney, Dehnen, & Bertelli [58] and Liu & Chaboyer [59]. In particular, 
we assume that the thin disk started forming 11.2 Gyr ago (which correspond 
to z ~ 2.1 in our cosmological model) and kept forming stars at the same 
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Fig. 6. Aitoff projection of the intensity of the Galactic neutrino background on the 
sky. Color coding scales logarithmically with the background intensity, and Galactic 
coordinates are used. Left and right panels refer to model Gl and G2, respectively. 

pace till to the present epoch. This is consistent with present data and with 
the assumption of a Salpeter IMF [58]. Note, however, that the slope of the 
initial mass function near 1 M & proves to be degenerate with the rate at which 
the star formation rate declines. We assume a synchronous formation of the 
galactic halo which dates back to 12.2 Gyr ago {z ~ 3), and a constant SFR 
in the thick disk lasting from 12.2 to 11.2 Gyr. As regards the bulge, recent 
HST/WFPC2 [60] and near-IR data [61] seem to suggest that it is coeval with 
the halo and that no significative amount of new stars formed after the intense 
initial starburst activity. 



Sky-maps of the Galactic neutrino background intensity are shown in Figure 
6. To facilitate a quantitative reading of the plots, in Figure 7, the back- 
ground intensity is plotted as a function of Galactic longitude for a number 
of ID-cuts at fixed latitude. 7 Both figures have been obtained by integrat- 
ing 1^ (r, E) over the particle energy, keeping r fixed. Note the spike in the 
background intensity towards the galactic centre, caused by the extreme stel- 
lar densities in the bulge. Independently of the Galaxy model, the amplitude 
of the neutrino background varies by roughly 3 orders of magnitude between 
the Galactic centre and the high-latitude isotropic component. Because of this 
large variation, all the different structural component of the Galaxy are clearly 
discernible from the maps. The amplitude of the high-latitude Galactic signal 
is always larger than the cosmic background (by a factor of a few in Gl and 
by roughly an order of magnitude in G2). This makes the detection of the 
cosmic background extremely difficult. 

Integrating the background over the whole sky, one obtains a neutrino flux of 
943 cm" 2 s" 1 , for model 1, and 333 cm" 2 s" 1 , for model 2. The corresponding 

7 Note that our maps have infinite angular resolution and should be convolved with 
an instrumental window function in realistic conditions. 
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Fig. 7. The Galactic neutrino background intensity as a function of direction of 
observation on the sky. Different lines show the t dependence at fixed latitude b. 
From top to bottom, b = 0°, 6°, 15°, 30°, 45°, 90°. Continuous lines refer to model 
Gl, while dashed lines are for G2. 
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Fig. 8. Energy spectrum of the Galactic neutrino background averaged over the 
whole sky. 

energy spectra are shown in Figure 8. In Table 3 we list the contributions 
to the flux from the different structural components of the Galaxy. At every 
latitude, the counts are dominated by neutrinos emitted in the Galactic disk 
which, independently of the Galaxy model considered, account for more than 
95 per cent of the total. 
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Table 3. Galactic neutrino background. 



Galaxy model 


Total 


Thin disk 


Thick disk 


Halo 


Bulge 




(cm- 2 s-') 


(cm- 2 !" 1 ) 


(cm-^" 1 ) 


(cm^s- 1 ) 


(cm-^" 1 ) 


1 


943 


923 


17.4 


1.72 


0.66 


2 


333 


314 


15.3 


1.12 


3.17 



4 Concluding Remarks 

We have estimated two neutrino backgrounds on the sky, in the 50 keV-10 
MeV window: the Galactic component from stars presently burning Hydrogen 
and Helium in the Milky Way and the cosmic contribution from thermonuclear 
reactions in other galaxies. On the grounds of the recent findings of Super- 
Kamiokande, SNO, and KamLAND, these neutrinos, born as v e , propagate as 
a superposition of (at least) two mass eigenstates. When (and if) interacting 
with matter, they will show up as neutrinos of all flavors, with comparable 
probabilities. 

Depending only slightly on some model details (namely, the assumed func- 
tional form of the SFR and the IMF), the background intensity of the cosmic 
component peaks at E ~ 0.1—0.3 MeV, where 1^ = 3— 7cm~ 2 s^ 1 MeV -1 sr" 1 , 
and sharply drops by 2 orders of magnitude above 1-2 MeV. The correspond- 
ing flux at the Earth's surface (integrated over all energies) ranges between 20 
and 65 cm" 2 s" 1 , while the associated energy density is ~ 5 x 10~ 4 eVcm -3 . 
Assuming that the sum of the masses of the three neutrino flavors is at least 
1/30 eV (as suggested by the oscillation patterns of solar and atmospheric 
neutrinos) implies that this energy density is negligibly small with respect to 
the contribution from the thermal (T Te \ = 1.95 K) background of relic neutri- 
nos originated in the early universe (w re i > 3.7eVcm~ 3 ). Similarly, the num- 
ber and entropy densities of cosmic thermonuclear neutrinos (n ~ 10~ 9 cm~ 3 , 
s ~ 10~ 7 /cbcm -3 ) are overwhelmed by the corresponding quantities for the 
primordial background (n re i — 10 2 cm -3 , s rc \ ~ 10 3 /cbcm -3 ). 

Thermonuclear neutrinos reaching the Earth today have been typically pro- 
duced by stars around z ~ 2. Even though they have been around for a 
significant fraction of the life of the universe, they hardly interacted with 
anything. As far as we know, the thermonuclear neutrino background is abso- 
lutely transparent to any diffuse radiation impinging onto it since the mean 
free path A = (na)' 1 exceeds the Hubble length for any reasonable value for 
the interaction cross section a. 

Attributing a fraction / of the observed 7-ray background in the sub-MeV re- 
gion to neutrino radiative decay, one can derive a lower bound for the radiative 
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lifetime, valid for any neutrino type and mass m, r^jm > (2 — 6) / 1 x 10 12 s 
eV" 1 . This can be easily translated into a limit for the electromagnetic form 
factors, which, after accounting for the contribution of active galactic nuclei 
(i.e. by assuming / ~ 1CT 2 ), is competitive with the results obtained from the 
analysis of the SN1987A neutrino burst. 

Independently of the direction of observation on the sky, the Galactic ther- 
monuclear background dominates its cosmic counterpart. Even at high lati- 
tudes, Galactic neutrinos outnumber the cosmological background by a factor 
which ranges between a few and 10 depending on the details of the Galaxy 
model (see Figures 6 and 7). The integrated flux of Galactic neutrinos over 
the whole sky ranges between 300 and 1000 cm~ 2 s _1 . The emission from stars 
in the Galactic disk contributes more than 95 per cent of the signal. 

Is there any prospect for detecting the Galactic component and, in the long 
term, to derive an observational neutrino map of the Galaxy? This clearly 
requires directional detectors, since one has first to distinguish the Galactic 
neutrino flux from the solar component. In the last ten years, large Cherenkov 
detectors have been built and successfully operated (Kamiokande, Super- 
Kamiokande and SNO), recording the directional signal of target electrons 
scattered from the impinging neutrinos. Analyzing the angular distribution of 
the events with respect to the Sun direction, the peak due to solar neutrinos 
is clearly discernible over a flat background [62]. So far, this method has been 
applied to the detection of the most energetic solar neutrinos originating from 
Boron decay (corresponding to a measured flux of 2.4 x 10 6 cm -2 s _1 ). In this 
case, the decay of 222 Rn in water, natural radioactivity, reactor antineutrinos, 
and radioactive decay of muon induced spallation products account for most 
of the instrumental background. After applying all possible cuts, background 
counts above 5 MeV are about ten times the number of events from solar neu- 
trinos. Even though this is not the optimal energy window for the detection 
of Galactic neutrinos, it is interesting to note that, for E u > 5 MeV, their flux 
amounts to ~ 1 cm~ 2 s _1 , which is nearly 7 orders of magnitude smaller than 
the instrumental background. This shows how far we are from their detection. 

Dedicated experimental designs and techniques should then be devised to 
detect cosmic and Galactic neutrinos. Let us consider, for instance, the cos- 
mic (anti)neutrino background from Type II supernovae. Kaplinghat et al. [5] 
found an upper bound to the flux of 54 cm -2 s _1 while, more recently, Ando, 
Sato, & Totani [6] estimate c-n ~ 11 cm -2 s _1 . Even though there is no energy 
region in which the contribution from supernova neutrinos dominates over the 
backgrounds, their presence might be detectable with Cherenkov telescopes as 
a distortion in the spectrum of electron and positrons produced in the decay 
of invisible muons. Ten years of integration at SuperKamiokande (or one at 
the proposed HyperKamiokande) should already guarantee a 1 a detection of 
the supernova neutrino background at ~ 20 MeV [6] . 
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Over the next decade, the focus of solar neutrino experiments will shift grad- 
ually from the relatively high energy 8 B neutrinos to the low energy, 7 Be, p-p, 
and CNO neutrinos. This is an encouraging perspective for the detection of 
a Galactic signal since, for E < 15 MeV, the number density of thermonu- 
clear neutrinos in the cosmic and Galactic backgrounds is expected to exceed 
that generated by Type II supernovae. Regrettably, instrumental backgrounds 
tend to increase in the MeV region. Therefore, in order to use the solar neu- 
trino experiments to measure a Galactic signal, it is mandatory to suppress 
the noise in the detectors by several orders of magnitude by using more effi- 
cient screening procedures and accurate selection of radiopure materials. Note 
that extremely long times of integration might be anyway necessary to detect 
Galactic neutrinos with high statistical significance. 

In summary, thermonuclear reactions in the Galaxy produce an anisotropic 
neutrino background in the MeV region. At the Earth's surface, this corre- 
sponds to an integrated flux of ~ 1000 cm~ 2 s -1 . Even though current detec- 
tors are unable to distinguish this signal from other backgrounds, it is not 
inconceivable that Galactic neutrinos might be revealed in the relatively near 
future with dedicated experimental setups. 
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A Entropy of a fermionic system 

Fermions have antisymmetric wavefunctions under particle exchange and this 
implies that no two fermions can exist in the same quantum state (Pauli 
exclusion principle). The number of ways an energy level Ei corresponding to 
gi degenerate quantum states can be populated with rii fermions is given by 
the binomial distribution with boxes, where rij contain 1 particle and gi — rij 
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are empty, i.e. 



9l - (A.l) 



nj. (ft -Ui)\ ' 

Allowing for a system with different energy levels, the total number of config- 
urations is 

We define the entropy of the system as 



S = k h lnW({n t })~ 
-hJ2{9i( ln 9i - 1) -rii(hini - 1) 



(A.3) 



(9i -nj)[ln(ft -m)-l]} , 



where, assuming both rii and g^ are ^> 1, we used Stirling's approximation 
Inn! ~ n(\nn — 1). The problem with eq.(A.4) is that it requires knowledge 
of the whole set {rii}, which, in practice, is a formidable task. A standard 
approximation is to replace W^({nj}) with W^({nj}) (where n« denotes the 
mean occupation number of the energy level Ej) in the above definition of 
entropy. In this case, neglecting 1 compared to gi, one gets 



rii In — 



(£-0--H)]- (A4 » 



It is convenient to introduce the mean occupation number of a single quantum 
state fi = fii/gi to obtain 




In (1 - fi) 



(A.5) 



If the energy levels of the system are closely spaced, one may replace the 
discrete probability density fi by a continuous function f(E). This is easily 
done for a set of free particles moving in an (essentially) infinite volume. In this 
case, all the sums over the energy levels J2i9i ■ ■ ■ rnay be replaced by integrals 
over the accessible phase-space volume f...dT/h 3 with dT = d 3 qd 3 p. For 
ultra-relativistic particles, for which E = cp, the entropy per unit volume is 
then given by: 



s = 



4?r k h g i 
c 3 h 3 



r°° ( 1 

I E2 i f{E Hm 



In [1 - /(£)] dE , (A.6) 



with g s the number of helicity degrees of freedom. This can be expressed in 
terms of the particle number density, 



n 



/ E 2 f(E) dE , 
J o 



as 



47T & 
c 3 h 3 

s = x k h n , 



(A.7) 
(A.8) 
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with 



_ g £ 2 {/(£) In [f(E)-i - 1] - In [1 - f(E)]}dE 
X J °°E*f(E)dE 1 • ) 



a dimensionless coefficient. 



B Photon Spectrum from Radiative Decays 



Let us assume that a set of cosmological sources produce neutrinos with a 
rate h v (z) (particles per unit co moving volume per unit time) and with a 
redshift dependent energy spectrum S(E v ,z) (normalized to unity). Then, 
the differential comoving neutrino number density per unit energy at redshift 
z is given by 



dn v 
dK 



(E„,z)= / h v {z')S 

J z 



'E v {l + z') \l + z' 



l + z 



l + z 



dt 



dz' 



H (1 + z) Jz 



S 



'E v {l + z') 
l + z 



dz' 
n v {z') dz 1 



(B.l) 



[n {i + z'f + n k {i + z'f + n K Yi 2 ' 



If neutrinos of energy E v decay radiatively with a (rest-frame) lifetime r 7 
Hq 1 , the decay rate in the cosmic frame will be n„(z) h(E u ) r 7 , where the 
Lorentz factor "({Ej) = E u /m u c 2 . Assuming that the secondary neutrino is 
very light compared to the parent one implies that the decay photons have 
E ~ E u /2, so that the decay processes will produce a photon background with 
present-day intensity: 



4tt Jo 



dz(l + z) 



00 dE v dn u ^ 
o l{E v ) r 7 dE v 



E v ,z) 5 D [E(l + z)- 



E v 



dt 



dz 



(B.2) 

Inserting equation (B.l) into equation (B.2) and integrating over the Dirac- 
delta distribution, one eventually gets 



■(F) - mv ° 2 1 f°° dz r°° ,h u (z')S(2E(l + z'),z') 

3{ } ~ 4^ HfE Jo (l + z) 2 K(l + z) J z K(l + z>) ' 

(B.3) 

with K(x) = [VLqx^ + VL k x 2 + Ha] 1 ^ 2 . In Section 2.7, we used this expression 
to constrain the neutrino decay time. 



Considering, for simplicity, a redshift-independent, monochromatic neutrino 
spectrum at emission, S(E V , z) = 5d(E u —E ), the neutrino energy distribution 
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reduces to 
dn 



dE v 



El'* h u (\E,(l + z) / E u ] - I) 



for E v < E 



H Q [fi El (1 + zf + Q k E v El (1 + zf + tt A £3]i/2 
for > £ 

(B.4) 

where cosmological expansion smeared the energy of neutrinos into a contin- 
uous spectrum. Similarly, for the photon background intensity, one gets 



c 1 m v c 



n„{~z{E)) 



for E v < E /2 



4tt 2 r 7 # 2 £ 2 K(E /2E) Ji w 2 K{w) 
for E v > E /2 

(B.5) 

with z(E) = (E /2E) — 1. For an Einstein-de Sitter universe (fi = 1, = 0), 
this becomes 



2 3/2 



4vr 5 




for E > E /2 
(B.6) 

with j(E) oc 22 x / 2 for E <^ Eq. 8 On the other hand, an empty universe 
(Qo = = 0) would give 



c 1 m^c 2 n u (z(E)) ( E t 



4tt 2 r 7 22 2 Eg V £ 




1 - 



2E\- 
~E~J 



for E < E /2 
for E > Eq/2 



(B.7) 



Finally, the effect of a non-vanishing cosmological constant can be sketched 
by looking at a de Sitter universe (fi = 0, ft a — 1) 



c l^nME)) ^* Cl-C^V) for £ < B./2 



4tt 2 r 7 




\E J 



for £ > E /2 . 

(B.8) 
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